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Geospatial relationships of air pollution and acute asthma 
events across the Detroit- Windsor international border: 
Study design and preliminary results 

Lawrence D. Lemke 1 , Lois E. Lamerato 2 , Xiaohong Xu 3 , Jason C. Booza 4 , John J. Reiners Jr 5 , Delbert M. Raymond III 6 , 
Paul J. Villeneuve 7,8 , Eric Lavigne 9 , Dana Larkin 2 and Helene J. Krouse 10 

The Geospatial Determinants of Health Outcomes Consortium (GeoDHOC) study investigated ambient air quality across the 
international border between Detroit, Michigan, USA and Windsor, Ontario, Canada and its association with acute asthma events in 
5- to 89-year-old residents of these cities. N0 2 , S0 2 , and volatile organic compounds (VOCs) were measured at 100 sites, and 
particulate matter (PM) and polycyclic aromatic hydrocarbons (PAHs) at 50 sites during two 2-week sampling periods in 2008 and 
2009. Acute asthma event rates across neighborhoods in each city were calculated using emergency room visits and 
hospitalizations and standardized to the overall age and gender distribution of the population in the two cities combined. Results 
demonstrate that intra-urban air quality variations are related to adverse respiratory events in both cities. Annual 2008 asthma rates 
exhibited statistically significant positive correlations with total VOCs and total benzene, toluene, ethylbenzene and xylene (BTEX) at 
5-digit zip code scale spatial resolution in Detroit. In Windsor, N0 2 , VOCs, and PM 10 concentrations correlated positively with 2008 
asthma rates at a similar 3-digit postal forward sortation area scale. The study is limited by its coarse temporal resolution 
(comparing relatively short term air quality measurements to annual asthma health data) and interpretation of findings is 
complicated by contrasts in population demographics and health-care delivery systems in Detroit and Windsor. 
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INTRODUCTION 

A growing body of epidemiological research has shown that 
ambient air pollution, and in particular traffic-related air pollution, is 
a contributing cause of premature mortality and morbidity. 1,2 Recent 
studies have reported adverse respiratory events from exposure to 
traffic emissions with intensified effects in children. 3-5 A systematic 
review of the research on long-term exposure to ambient air pollu- 
tion and chronic diseases also suggests an increased morbidity and 
mortality in adults from lung cancer, cardiovascular diseases, and 
respiratory conditions. 6 Major urban air pollutants such as nitrogen 
dioxide (N0 2 ), a recognized marker of traffic-related pollution, 7 have 
been associated with higher rates of respiratory morbidity, 8 particu- 
larly asthma. 9,10 Few studies, however, have evaluated the spatial and 
temporal variability of air quality and its effects on health conditions 
between adjacent international cities sharing a common airshed. The 
international context of this study provides unique opportunities to 
compare population health outcomes in response to environmental 
conditions in adjoining cities with different environmental regula- 
tions, public health policies, and health-care delivery systems. 

Epidemiological analyses of the health-related effects of air 
pollution are complicated by the need for accurate characterization 
of the effects from exposure to airborne pollutants (e.g., increased 



asthma exacerbations) and the availability of information indicative 
of exposure (e.g., ambient air quality measurements). Accurate 
characterization of exposure gradients for all members of a given 
study population is particularly challenging because routine 
monitoring networks rarely have enough sites to capture small- 
scale spatial variability within a city. The significance of monitoring 
resolution is emphasized by findings that urban N0 2 concentrations 
measured 50 m apart can vary by as much as 300% n and that 
traffic-related pollutants fall off to background levels within 300 m 
of major roadways. 12,13 Additionally, recent studies suggest that 
intra-urban exposure to airborne pollutants may be associated with 
larger health effects than inter-urban variability. 14,15 Consequently, 
there is a need to characterize intra-urban gradients for pollutants 
at neighborhood scales particularly for studies focused on 
characterizing adverse health effects associated with long-term 
exposure. This can be accomplished by establishing short-term, 
multi-site measurement networks (i.e., saturation monitoring 16 ). 

In addition, air pollutant concentrations also vary over time so that 
exposure is a function of both 'where' you are as well as 'when' you 
are there. Temporal variability has been documented on hourly, 
diurnal, and seasonal time scales, each of which may be important 
for asthma exacerbation or prevalence. For example, significant 
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associations have been documented between high pollutant 
concentrations and asthma hospitalizations 17 and other respira- 
tory symptoms 18 with time lags ranging from 0 to 5 days. In another 
analysis, temporal variability exceeded spatial variability for widely 
spaced measurements of 39 volatile organic compounds, carbonyl 
compounds, trace elements, and extractable organic mass species 
in an urban area. 19 Consequently, the questions of how, when, and 
where to measure or estimate exposure to airborne contaminants 
are of fundamental importance when assessing health outcomes. 

Air Quality Studies in Detroit and Windsor 
Detroit, Michigan and Windsor, Ontario, are neighboring cities 
separated by the Detroit River along the USA-Canada border. The 
economies of both cities share strong ties to the automotive 
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industry along with high air pollution levels relative to most North 
American cities. Concerns about air quality have motivated recent, 
separate exposure studies in each city. In Windsor, land-use 
regression (LUR) models were developed for N0 2 , S0 2 , benzene, 
and toluene concentrations measured with a dense network of 50 
monitors during 2-week sampling periods in each of the four 
seasons of 20 04. 20,21 Later, the Windsor, Ontario Exposure 
Assessment Study (WOEAS) utilized personal, indoor and outdoor 
monitors to estimate the exposure of adults and asthmatic children 
living in Windsor to air pollutants in 2005 and 2006. 22 Most 
recently, positive associations were found between emergency 
department visits for asthma and air pollution measured at central 
monitoring stations in Windsor from 2002 to 2009 23 

In Detroit, the Detroit Air Toxics Pilot Study demonstrated 
variability of particulate matter and air toxics measured at eight 




Figure 1. Detroit/Windsor study area with sampler locations and postal code areas. 



Table 1. 


Instrumentation, sample preparation and laboratory analysis. 






Analyte 


Sampling 
type 


Sampler 


Sampler preparation 


Laboratory analysis 


N0 2 and S0 2 Passive 


Ogawa diffusion samplers 


Pretreated filter obtained 


Environment Canada (Egbert, Ontario) 








from Ogawa 


using ion chromatography 


VOC 


Passive 


3M #3500 organic vapor monitor 


Obtained from 3M 


Airzone One (Mississauga, Ontario), gas 
chromatography/mass selective 
detector with species dependent 
detection limits, 26 VOCs 


PAH 


Active 


URG Model URG-2000-25AB Personal 


PUF samplers cleaned and 


Airzone One (Mississauga, Ontario), gas 






Pesticide Samplers; 


pre-conditioned, Airzone 


chromatography/mass selective 






SKC Leland Legacy pump at 2.0 Lpm 


One 


detector with species dependent 






flow rate 




detection limits, 16 PAHs and 7 
pesticides 


PM 


Active 


3-stage Harvard Cascade Impactors; 


Teflon filter and PUFs 


Alberta Research Council, gravimetric 






SKC Leland Legacy pump at 5.0 Lpm 


conditioned before and after 


method with lab detection limit of 






flow rate; 


sampling, Alberta Research 


4(ig/filter 






PMt on a Teflon filter, PM 1 _ 2 . 5 and 


Council (Vegreville, Alberta) 








PM25-10 ° n individual Polyurethane 










Foams (PUFs) 
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Table 2. Sampler retention and species included in September 2008 and June 2009 analyses. 



Pollutant 



September 5-20 2008 



May 29-June 13 2009 





Samplers 
deployed 


Samplers 
retrieved 


Samples 
accepted 


Species included in 
analysis 


Samplers 
deployed 


Samplers 
retrieved 


Samples 
accepted 


Species included in 
analysis 


N0 2 


100 


99 


98 


Included 


100 


100 


93 


Included 


S0 2 


100 


99 


98 


Excluded 


100 


100 


90 


Included 


VOCs 


100 


100 


98 


1 2 out of 26 


100 


100 


99 


9 out of 26 


PAHs 


50 


50 


38 


16 out of 23 


50 


50 


40 


8 out of 23 


PM 


50 


50 


38 


Included 


50 


50 


38 


Included 



Table 3. Detroit and Windsor population demographics (%). 





Detroit 
201 Cf 


HFHS Detroit 
2008 b 


Windsor 
2006 c 


Male 


47.3 


39.2 


48.8 


Female 


52.7 


60.8 


51.2 


Age 5-14 


14.3 


10.9 


12.3 


Age 15-44 


41.9 


48.5 


43.3 


Age 45-64 


25.4 


26.3 


24.1 


Age 65 + 


11.4 


7.9 


14.3 


Total 


93.0 


93.7 


94.0 


White d /not a visible 


10.6 


12.3 


79.0 


minority 6 








Black or African 


82.7 


78.6 


3.9 


American d /Black e 








Other 


6.7 


9.2 


17.1 



a Source: US Census 2010. 

b HFHS data exclude patients younger than 5 or older than 89. 
c Source: Statistics Canada 2006. 
d US Census classification 
Statistics Canada classification. 



locations in the Detroit metropolitan area in 2001 and 2002. 
Subsequently, the Detroit Exposure and Aerosol Research Study 
(DEARS) compared air pollutant concentrations measured at a 
central monitoring site with indoor, outdoor, and personal moni- 
toring measurements distributed throughout six Detroit neighbor- 
hoods from 2004 to 2007. 25 Although the spatial resolution of the 
outdoor measurements was limited, the DEARS study found 
spatial variability and gradients in ambient concentrations of 
particulate matter 26,27 and volatile organic compounds 28 Modest 
associations were found between community-based and personal 
exposure measurements of N0 2 in the DEARS study. 29 

Study Objectives 

The Geospatial Determinants of Health Outcomes Consortium 
(GeoDHOC) study was designed to explore associations between 
short-term ambient air pollution concentrations and acute asthma 
events at a high level of spatial resolution. Our approach is 
motivated by the expectation that correlations among mappable 
environmental attributes and health indicators can be used to 
better understand and manage urban community health. Specific 
aims included (1) collection and modeling of simultaneous air 
pollutant concentrations in Detroit and Windsor, (2) collection of 
concurrent data on asthma prevalence using asthma-related 
patient encounters, and asthma exacerbations using asthma- 
related emergency department visits and hospital discharge 
records in Detroit and Windsor, (3) integration of the environ- 
mental and asthma data into a geographic information system 
(GIS) database, and (4) analysis of spatial associations among air 
quality and asthma metrics. This paper provides an overview of the 
GeoDHOC project by describing the study design, implementation 
strategies, and initial findings for spatial models relating pollutant 
distributions to acute asthma events in Detroit and Windsor. 



METHODS 

Geospatial Database 

A geospatial database was developed to manage the large amount of 
information needed for this study. Data assembled for use in sampler 
allocation and land-use regression modeling included population character- 
istics, land-use characteristics, proximity-to-point sources in the US 
Toxics Release Inventory 30 and the Canadian National Pollutant Release 
Inventory, 31 topography, traffic volumes, and road networks. These public 
data were obtained through local, state/provincial, and national sources 
including the Southeast Michigan Council of Governments, US 
Environmental Protection Agency, US National Oceanic and Atmospheric 
Administration, City of Windsor GIS Division and Infrastructure Services, 
Essex Region Conservation Authority, Statistics Canada, and Environment 
Canada. Together with air sampling and health data, these data sets were 
integrated into a GIS to facilitate data visualization and geospatial modeling. 

Data rectification. Configuration of disparate data sets across the 
international study area (Figure 1) posed a significant database manage- 
ment challenge. Despite the proximity of Detroit and Windsor, many of the 
geographic data files used in this study had different geospatial 
projections (UTM versus State Plane), measurement systems (imperial 
versus metric), data collection methodologies (e.g., traffic counts, emissions 
inventories), scalar resolutions, or typologies (e.g., land use). In addition, 
data files were formatted for use in different software programs (e.g., 
Maplnfo versus ESRI ArcGIS). 

If data from one city were more complete or of a higher resolution, they 
were reduced to the standards of the lesser quality data for comparative 
purposes. In many cases, data quality was high for both cities, but data sets 
contained different typologies. For example, road and land-use classifica- 
tions are different in Detroit and Windsor. Thus, data categories for both 
cities were reclassified to create a new merged standard without 
deterioration in data quality. In other cases, such as point source emissions 
data, differences in reporting conventions between the US and Canada 
required us to keep the information in its original format. 

Air Quality Sampling Methods 

Air sampling involved a combination of 50 paired active/passive and 50 
passive-only devices deployed at 100 sampling sites. 32 Two sampling 
rounds of 2-week duration were conducted, the first in September 2008 
and the second in June 2009. Measured pollutants included N0 2 , S0 2 , PMt 
(particulate matter <1/^m in aerodynamic diameter), PM q _ 2 . 5 , PM 2 . 5 -io, 
and selected volatile organic compounds (VOCs) and polycyclic aromatic 
hydrocarbons (PAHs). These pollutants were selected for reasons including 
their known influence on respiratory health, their relationship to land use 
and traffic, and investigator experience measuring them in Windsor and 
other cities. 

Air sampling equipment. Air sampler specifications are given in Table 1. 
Equipment included 50 sets of active air monitoring equipment using SKC 
Leland Legacy sampling pumps (SKC Inc., Eighty Four, PA, USA) powered 
by an AC adapter and housed in weatherproof cases. The pumps ensured 
sufficient air flow to operate Harvard Cascade Impactor particulate 
monitors (Harvard School of Public Health, Boston, MA, USA) and URG 
personal pesticide samplers (model URG-2000-25AB, Chapel Hill, NC, USA). 
The latter were used to measure PAHs. Additional sampling equipment 
included 100 sets of passive diffusion monitors consisting of a rain shield 
and mounting bracket to hold the Ogawa passive badges (Ogawa & Co, 
Pompano Beach, FL, USA) used to monitor N0 2 and S0 2 , and the 3M Model 
3500 Organic Vapor Samplers (3M Canada, London, Ontario) used to 
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measure 26 different VOCs, including benzene, toluene, ethylbenzene, and 
xylene (BTEX). 

Sampler allocation. Existing data on transportation routes, population 
density, and the location of suitable sampling facilities were analyzed to 
design a sampler allocation plan for the study. Constraints included the 
number of available air sampling units and the need for safe and secure 
locations with a continuous power supply for active samplers. Methods 
employed by others including location-allocation modeling, 33 criteria- 
based targeted sampling, 34 and uniform distribution sampling 35 were 
considered. However, the absence of requisite mapped pollution gradients 
for Detroit and an at-risk population gradient for Windsor precluded 
location-allocation and criteria-based targeted sampling. Consequently, a 
site-specific methodology was developed for this study. 

Initially, we used land area and city shapes to determine that a 
DetroitWindsor sampler ratio of 34:1 6 would optimize the desired spatially 
stratified sampling design. Allocation of individual samplers then 
proceeded in two steps to accommodate the need to locate paired 
active/passive and passive-only samplers. Based on our need for a 
continuous power supply and security concerns, active/passive samplers 
were placed at fire and police stations in Detroit and municipal and 
nonprofit institutions in Windsor. Specific sites were selected to lie as close 
as possible to the centroid of a 3 x 3 km cell grid pattern overlain on each 
city. Passive samplers were then allocated within the same 3 km grid, with 
first priority given to cells that did not contain a paired active/passive 
sampler. Secondary priority was given to maximizing distance from other 
samplers, minimizing proximity to schools, and collocation of samplers 
with permanent governmental air monitoring sites in Detroit and Windsor 
and six sites tied to prior studies in Windsor. 21 Collocated sites were used 
for measurement validation. 32 The average spatial resolution of the 
sampling design across the study area is 1 measurement per 5 km 2 for 
passive samplers, and 1 measurement per 10 km 2 for active samplers. 

Preparation and transport of air samplers. Measures were undertaken to 
ensure consistency and reliability in sampling across the international 
airshed. Monitoring equipment and samplers were prepared at the 
University of Windsor following rigorous quality control procedures. 
Samplers with prepared filters were transported in sealed bags inside 
coolers to maintain a constant temperature and prevent contamination or 
off gassing in transit prior to deployment at predetermined sampling sites 
in Detroit and Windsor. Detroit samplers were carried across the 
international border by private automobile (as opposed to commercial 
transport) so that study investigators could remain with the samples to 
avert contamination during border inspection. Customs officials on each 
side of the border were notified in advance of shipments and X-ray 
inspection was requested to avoid exposing filters to ambient air and 
temperature changes (both prior to deployment and following air sample 
collection). Transportation by car also helped minimize delays inherent in 
commercial international shipping. 

Sampler deployment and retrieval. Samplers were distributed between 
Detroit and Windsor (Table 2) based on the allocation method described 
above. Combined active and passive samplers were deployed at 50 
locations on the first day of each sampling round, and 50 additional 
passive-only samplers were deployed on the second day. Samplers were 
deployed on light poles, downspouts, and tripods at predetermined 
specified heights above the ground in both cities. 32 All samplers were 
retrieved 2 weeks after their initial deployment. Municipal and logistical 
constraints necessitated the movement of six sampling sites between the 
September 2008 and June 2009 sampling events. Relocated sites were 
placed within 200 m of the original locations. Four additional passive sites 
were added in the 48217 zip code area in 2009 to extend coverage within 
southwest Detroit. 

Quality assurance. Quality assurance (QA) and quality control (QC) 
measures were applied to ensure consistency, reliability, and reproduci- 
bility of the study results. Field equipment was tested and calibrated prior 
to deployment. Samples were handled and stored by trained technicians 
using written standardized operating procedures to enforce sampling 
consistency and minimize potential contamination. Field and laboratory 
blanks were deployed to quantify the sample mass attributable to 
handling and transportation. Because of the many samples and sampling 
stages involved, an internal sample manifest, nomenclature, and tracking 
procedure was implemented including log sheets and chain of custody 
records. Deviations from the established protocol were noted and log 
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sheets were reviewed for discrepancies or unexplained observations. After 
retrieval, samplers with out-of-tolerance flow rates, insufficient pump run 
times, or evidence of tampering were excluded from further analysis 
(Table 2). In addition, individual species with <80% of sites registering 
values above laboratory detection limits were excluded from further 
analysis (Table 2). A more detailed description of sampling methods, 
precision of instruments, and QA/QC procedures employed is given by 
Miller et al. 32 and Wheeler et al. 22 



Health Data Census Methods 

Study population. Individuals 5-89 years of age living in the Detroit and 
Windsor metropolitan areas are the population of interest for this study. 
Asthma is the health outcome of interest. Children under the age of 5 were 
excluded because objective measures of lung function such as spirometry 
are standardized only for individuals 5 years of age and older, and it is 
often difficult to differentiate asthma from other respiratory conditions in 
very young children. 36 

Patient demographics. As of 201 0, Detroit had ~ 71 4,000 residents; ~ 30% 
were under 20 years of age and 11% were 65 years and older. 37 In 2008, the 
Henry Ford Health System (HFHS) provided health-care services to 
> 125,000 Detroit residents. The demographics of HFHS patients living in 
Detroit generally reflect those of Detroit— the majority of the HFHS Detroit 
patient population is African American, but a larger proportion of HFHS 
patients are female (Table 3). HFHS patients residing in Detroit had ~400 
hospital admissions and 2,500 emergency room visits for asthma in 2008. 



Windsor, similar to Detroit, is a diverse city that is home to many ethnic 
and racial groups. According to the 2006 census, 38 ~25% of Windsor 
residents are under 20 years of age and 14% of the population is 65 years 
and older (Table 3). In contrast to Detroit, where ~80% of the population is 
African American, a small minority (4%) of Windsor residents are black and 
other ethnic and racial groups comprise 17% of the population (Table 3). 
Windsor Regional Hospital (Metropolitan Campus and Western Campus) 
and Hotel Dieu Grace Hospital are the main medical facilities. Together, 
their campuses experienced 423,000 ambulatory visits in 2007/2008. 39 
Hospitalization data are captured in two national databases maintained by 
the Canadian Institute for Health Information (CIHI): the Discharge Abstract 
Database (DAD) 40 for hospital admissions, and the National Ambulatory 
Care Reporting System (NACRS) 41 for emergency room visits. In order to 
avoid duplicate counts (i.e., a patient admitted to the emergency room and 
subsequently admitted to the hospital), we used a unique patient 
identification number to allow tracking across databases and removed 
those that were captured in both databases for the same asthma event. As 
health care in Canada is available through a publicly funded system, the 
demographics of Windsor residents are assumed to be representative of 
patients receiving care in the Windsor medical facilities. 

Health outcome variables. Asthma measures collected as part of the 
greater GeoDHOC study include prevalence rates and acute asthma events. 
Asthma prevalence, defined as the number of affected individuals 
normalized to the serviced population residing within each environmental 
assessment area, was determined using diagnostic coding from clinical 
encounters. Acute asthma events were assessed by enumerating both 
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hospital admissions and emergency department visits with primary diag- 
noses of asthma. Service date and the appropriate ICD9-CM codes (493.xx) 
or ICD-10 code (J45) were used to identify health data for analysis. 
Descriptive analyses were used to characterize asthma morbidity in the 
population by age, gender, race, residence, and other sociodemographic 
variables. Geocoding by residential address was used to assign a spatial 
location to health outcome event records in Detroit, and Windsor records 
were assigned to the appropriate postal forward sortation area based on 
residential address. Acute asthma counts within each postal code area 
were stratified by age and gender to account for spatial differences in the 
underlying age-sex distribution of the study population. Institutional 
Review Board approval for the study was obtained from the participating 
institutions and all personal and health information was de-identified and 
coded to protect the identity of individuals. 



Spatial Correlation of Air Quality and Health Data 
Interpolated models of air pollutant concentrations across Detroit and 
Windsor were generated on a 300 x 300 m grid using ordinary kriging with 
variogram parameters derived for the entire Detroit-Windsor study area. 32 
Contour maps were created from these grids to depict the spatial 
distribution of measured analytes, including PM 2 . 5 and PM qo , which were 
obtained by summing grids for the measured PM q plus PM q _ 2 . 5 , and PM 2 . 5 



plus PM 2 .5_io components, respectively. Average concentration values for 
postal zones (5-digit zip code tabulation areas in Detroit and 3-digit 
forward sortation areas in Windsor) were calculated from the kriged models 
using spatial analyst tools in ArcGIS. The density of asthma occurrences, 
geocoded by patient residence address and normalized to patient 
population (i.e., persons with residential addresses in Detroit receiving 
health care from HFHS during 2008 and total population of Windsor in 
2006), was determined for individual postal zones in Detroit and Windsor. 
Pearson correlation coefficients and regression coefficients were calculated 
to assess the strength of associations between rates of acute asthma 
events and corresponding average concentrations of measured pollutants 
in each postal zone. Zip codes 48203 (Highland Park) and 48212 
(Hamtramck) were excluded in Detroit and forward sortation area N8N 
was excluded in Windsor because these postal zones include significant 
numbers of residents living outside the city boundaries. In addition, health 
information was not available for forward sortation area N8V. 



RESULTS 

Ambient Concentrations and Distributions 

Aggregated data for selected 2008 and 2009 air pollutant 

concentrations across the two metropolitan regions are presented 
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in Table 4. A more detailed description and analysis of the 
September 2008 measurements and spatial distributions of airborne 
contaminants was reported by Miller et al. 32,42 Mapped distributions 
of September 2008 N0 2 , total VOCs, total BTEX, PM 1# PM 2 . 5/ and 
PM 10 concentrations are shown on Figure 2. Spatial variability is 
readily apparent between and within both cities. Concentrations of 
the September 2008 mapped pollutants averaged across each 
postal zone using ArcGIS are shown in Figure 3. 

An independent sample f-test was used to evaluate the 
statistical significance of differences between mean concentra- 
tions measured in Detroit versus Windsor for each analyte. In 2008, 
the mean concentration of N0 2/ BTEX, VOCs, PAHs, PM V and PM 10 
were significantly higher (P<0.05) in Detroit than Windsor. In 
2009, mean concentrations of N0 2 , BTEX, and VOCs were higher in 
Detroit and the mean concentration of S0 2 was higher in Windsor 
(Table 4). 

Comparisons of concentration variance suggest greater varia- 
bility for some analytes within Windsor. The coefficients of 
variance (CV), obtained by dividing the standard deviation by 
the mean (Table 4), were consistently greater in Windsor for N0 2 , 
PAHs, PM 1# PM 25 , and PM 10 in 2008 and 2009, although CVs for 
VOCs were lower in Windsor. 

Health Outcomes and Correlations with Ambient Concentrations 
In Detroit, HFHS asthma-related health-care visits were appor- 
tioned by zip code tabulation area corresponding to each patient's 
residential address (Table 5). Approximately 3000 hospitalizations 
and ER visits were recorded at HFHS facilities for primary diagnosis 
of asthma for patients residing in the study area during 2008. In 
Windsor, asthma-related health-care visits quantified using 
hospital discharges and emergency room visits were apportioned 
to Canadian forward sortation areas corresponding to the 
patients' residential address (Table 6). The density of asthma 
occurrences for individual postal zones in Detroit and Windsor is 
shown in Figure 4. 

Visual comparison of Figure 4 with September 2008 air 
pollutant concentrations averaged by postal zone (Figure 3) 
suggests potential relationships between the standardized num- 
ber of asthma events and concentrations of some pollutants in 
Detroit and Windsor. With the exception of PM 25 and PM 10 in 
Detroit and PMt in Windsor, positive correlations are apparent in 
scatter plots of age-gender-standardized asthma events versus 
average pollutant concentrations by postal zone for Detroit 
(Figure 5) and Windsor (Figure 6). Although these plots display a 
fair degree of scatter, correlation and regression slope coefficients 
(Table 7) demonstrate moderate to strong postal code scale 
associations between standardized asthma events and N0 2/ VOCs 
and PM 10 in Windsor, as well as BTEX and VOCs in Detroit. 



DISCUSSION 

The observed association between N0 2 and acute asthma events 
in Windsor adds to a growing number of studies linking asthma 
morbidity to higher levels of N0 2 in other cities. Using coarse 
measures of N0 2 to estimate exposure, Hajat et al. 43 found an 
association between N0 2 and general practitioner consultations 
for asthma in London, England. More recently, Ko et al. 17 reported 
a significant association between asthma hospitalization rates and 
N0 2 in Hong Kong while Ho et al 44 found a significant relationship 
between adolescents' asthma attack rates and N0 2 in Hong 
Kong. A number of studies citing the production of N0 2 as a 
byproduct of roadway traffic also support an association between 
roadway and other point/non-point source proximity and 
asthma outcomes. For example, Oyana et al. 45 documented 
the clustering of adult asthma outpatient visits and hospitaliz- 
ations near roadways and a significant bridge complex in 
Buffalo, New York. Similarly, Newcomb and Li 46 concluded that 



Journal of Exposure Science and Environmental Epidemiology (2014), 346-357 



© 2014 Nature America, Inc. 



Detroit-Windsor air pollution and asthma 
Lemke et a I 



353 

Table 6. Calendar year 2008 asthma hospital discharges/emergency department visits and September 2008 average N0 2 , BTEX, PMt, PM 2 . 5/ and 
PM 10 levels for Windsor forward sortation areas. 



Forward 


Area 


Population 


2008 asthma 


2008 asthma 


Age-gender 


Sept 2008 


Sept 


Sept 


Sept 


Sept 


sortation 


(km 2 ) 


(2006 


hospitalizations 


exacerbations 


standardized 2008 


average 


2008 


2008 


2008 


2008 


area 




census) 


and ER visits 


per WOO pop a 


asthma 


N0 2 


average 


average 


average 


average 












exacerbations 


(p.p.b.) 


BTEX 

(fig/m 3 ) 


PM 1 

(fig/m 3 ) 


PM Z5 
(fig/m 3 ) 


PM 10 
(fig/m 3 ) 


N8P 


6.7 


9,270 


10 


1.08 


1.8 


9.7 


5.24 


5.2 


7.9 


10.1 


N8R 


4.6 


12,350 


22 


1.78 


3.2 


10.4 


4.85 


5.0 


7.8 


9.9 


N8S 


7.9 


22,845 


53 


2.32 


9.1 


11.5 


5.31 


5.2 


7.8 


10.8 


N8T 


9.9 


1 8,245 


50 


2.74 


7.9 


1 1.8 


5.78 


5.0 


6.8 


10.6 


N8W 


16.9 


24,810 


73 


2.94 


11.1 


13.3 


6.29 


5.0 


6.7 


10.9 


N8X 


9.6 


16,115 


77 


4.78 


12.8 


14.2 


7.65 


5.0 


7.0 


11.4 


N8Y 


7.5 


19,720 


78 


3.96 


13.5 


12.5 


6.21 


5.2 


7.3 


11.4 


N9A 


6.0 


26,140 


120 


4.59 


21.0 


15.0 


6.22 


5.4 


8.7 


14.2 


N9B 


7.6 


16,190 


68 


4.20 


10.6 


15.4 


8.46 


4.7 


6.7 


12.0 


N9C 


15.4 


12,495 


44 


3.52 


7.3 


13.9 


5.85 


5.0 


6.8 


13.6 


N9E 


13.6 


20,865 


25 


1.20 


3.1 


13.3 


6.92 


4.9 


6.5 


11.0 


N9G 


9.9 


1 7,505 


25 


1.43 


3.2 


13.0 


4.73 


5.6 


7.2 


12.0 


Windsor 


116 


259,240 


645 


2.49 




12.8 


5.9 


5.1 


7.3 


11.5 



denominator includes total population of Windsor in 2006. 




Figure 4. Year 2008 asthma events per 1000 population by postal code area. Note different scales for Detroit and Windsor. 



proximity to major roadways is a strong predictor of emergency 
department admission or hospitalization among children in 
Terrant County, Texas. In contrast, Lewis et al. 47 did not find an 
association between roadway proximity and physician-diagnosed 
asthma or asthma medication in east England, and Meng et al. 48 
reported the lack of a significant relationship between poorly 
controlled asthma in adults and N0 2 measurements from govern- 
ment air monitoring stations in Los Angeles and San Diego, 
California, USA. 

The observed associations between asthma and VOCs in Detroit 
and Windsor contribute to a small number of studies examining 
relationships between asthma and VOC concentrations, although 
comparisons to published studies are complicated by differences 
in how VOC data are reported. Weichenthal et al 49 reviewed six 
studies investigating indoor VOC exposure and childhood asthma. 
Three of those studies included only formaldehyde, which was not 



included in our analysis. The results for the other three studies 
were mixed with two reporting significant asthma risk associated 
with VOCs and one finding no increase in asthma risk or morbidity 
linked to total VOCs. Arif and Shah 50 found that physician- 
diagnosed asthma and the presence of wheezing in US adults was 
significantly related to increased exposure to separate groupings 
of four aromatic compounds (BTEX) and two chlorinated 
hydrocarbons (tetrachloroethene and trichloroethene) measured 
with passive personal exposure monitors. Kim et al. 51 studied 
biomarkers of VOC exposure (i.e., hippuric and muconic acid), and 
their relation to childhood asthma. They reported a significant 
increase of muconic acid in subjects with asthma. 

With respect to particulate matter and acute asthma events, our 
results are mixed. We found a significant association between 
asthma and PM 10 in Windsor, but correlations for other size 
fractions were not significant in either city. Although comparison 
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Sept 2008 Avg N02 (ppb) Sept 2008 Avg PM 1 (ng/m3) 




Sept 2008 Avg Total VOC (ug/m3) Sept 2008 Avg PM2.5 (p-g/m3) 




Sept 2008 Avg Total BTEX (ug/m3) Sept 2008 Avg PM10 (ug/m3) 



Figure 5. Detroit 2008 calendar year age-gender standardized HFHS asthma rate versus average September 2008 concentration by zip code 
for: (a) N0 2 ; (b) VOC; (c) BTEX; (d) PIVh; (e) PM 2 . 5 ; and (f) PM 10 . Lines show best fit linear regression. 



is sometimes confounded by inconsistent reporting of PM size 
fractions in the literature, results of other studies are also mixed. 
Meng et al. 48 demonstrated a relationship between poorly 
controlled asthma and PM 10 measured at government air 
monitoring stations in Los Angeles and San Diego, California, 
primarily in women; however, significant associations with PM 2 . 5 
were not observed. Hajat et al. 43 found that PM 10 concentrations 
measured at a central monitoring site were positively related to 
asthma treatments by general practitioners in London. Ko et al. 17 
also reported that PM 2 . 5 and PM 10 were both significantly 
associated with asthma hospitalization rates in Hong Kong. Con- 
versely, although trends in their data suggested a relationship 
with lung function, Tang et al. 52 reported no statistically significant 
results linking PM n , PM 25 , or PM 10 (assessed with personal 
monitors) to peak expiratory flow rate for asthmatic children in 
Sin-Chung City, Taiwan. Slaughter et al. 53 also reported a lack of 
statistically significant associations among four PM size fractions 
and hospitalizations or ER visits for asthma in Spokane, 
Washington. 

The GeoDHOC study has several important limitations. First, it 
relies upon ambient air pollutant concentrations for exposure 
estimates and does not address personal exposures of the study 
population within individual residences or places of occupation. 
Second, the collection of integrated air pollutant measurements 
over 2-week sampling periods misses daily or diurnal fluctuations 



in air pollutants such as traffic-derived N0 2 . This coarse temporal 
resolution was a necessary tradeoff of the study design 
incorporating a large number of samplers to increase spatial 
resolution. Third, the results reported here compare pollution 
measurements in September to asthma events recorded through- 
out a calendar year without accounting for seasonal variation. 
Prior studies in Windsor examining air pollutant variability across 
four seasons (summer, winter, spring, and fall) in 2004-06 found 
significant correlations among all seasons for N0 2 , S0 2 , and BTEX, 
suggesting a reasonable degree of consistency in source 
distributions for N0 2 and VOCs throughout the year. 21,54 Fourth, 
we drew acute asthma events from two different health databases 
representing populations from different countries and health-care 
delivery systems. The results reported here do not explicitly 
account for sociodemographic differences between Detroit and 
Windsor, nor do they account for potential differences in 
diagnoses and coded services between the populations. 
Moreover, by relying on hospital visit administrative data only, 
the study does not take into account the role that management of 
disease may have in exacerbation. Significant differences in any or 
all of these factors potentially could engender ecological fallacies 
arising from individual-level conclusions drawn from aggregate 
health effect estimates. 55 

One noteworthy difference between Detroit and Windsor is the 
number of acute asthma events observed per 1000 population. 
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Figure 6. Windsor 2008 calendar year age-gender standardized asthma rate versus average September 2008 concentration by forward 
sortation area for: (a) N0 2 ; (b) VOC; (c) BTEX; (d) PM 1# - (e) PM 2 . 5 ; and (f) PM 10 . Lines show best fit linear regression. 



Table 7. 


Linear regression slope (|31 values) and Pearson correlation coefficients (r) for air pollutant concentrations and asthma events by postal zone 


in 2008 (data in Tables 5 and 6). 
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The average population-normalized asthma events for HFHS 
patients in Detroit zip codes (Table 5, column 6) is ~10 times 
greater than the average number of asthma events in Windsor 
forward sortation areas (Table 6, column 5). Part of this 
discrepancy is attributable to the difference in denominator 
specification between the two cities (i.e., HFHS patients receiving 
health care in Detroit versus total population of Windsor). 
However, rates for Detroit are consistent with recent reports of 
5.1 asthma hospitalizations per 1000 Detroit residents, and 33.6 



asthma emergency department visits per 1000 Detroit children 
(<18 years) in 2004-2006. 56 Moreover, HFHS utilization rates for 
Detroit residents in 2008 are also consistent with rates for 2007 
and 2009, indicating that 2008 is not anomalous. Therefore, other 
factors inherent to the complexity of comparing health outcomes 
between cities in different countries are also likely to contribute to 
this disparity. For instance, differences in demographics including 
contrasts in income and ethnic and racial composition are readily 
apparent between Detroit and Windsor. In addition, Canada uses a 
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universal health-care model while the USA has a multiplicity 
of payer methods and a significant proportion of uninsured 
residents. HFHS is a private, nonprofit organization serving both 
commercially insured persons (largely employer-paid plans), 
government-financed insured persons (Medicare for the elderly 
and Medicaid for the indigent), as well as a significant portion of 
uninsured persons. Records indicate 19% of the HFHS asthma 
patients in this study received uncompensated care, and an 
additional 30% relied upon non-commercial insurance with 
variable benefits (e.g., uncovered ambulatory care), potentially 
elevating emergency care utilization rates in Detroit. 

Beta values for BTEX and VOCs in Detroit exceed those for N02, 
VOC, and PM10 in Windsor (Table 7), indicating a potentially 
greater increment of standardized acute asthma events per unit 
increase in ambient concentrations in Detroit. Conversely, higher 
correlation coefficients between standardized asthma events and 
pollutant concentrations were observed in Windsor. Nevertheless, 
the 95% confidence range of r values overlaps for all significant 
associations (Table 7), indicating r values are not distinguishable 
from each other or between the two cities. Hence, although 
average concentrations of air pollutants found in Detroit were 
generally higher than those found in Windsor (Table 4), the larger 
slope factors for BTEX and total VOCs observed in Detroit (Table 7) 
do not necessarily indicate stronger associations with asthma 
exacerbations in Detroit compared with Windsor. Rather, these 
initial results underscore the need for more rigorous investigation 
of how differences in population demographics, access to health 
care, and asthma management practices help to explain the 
observed differences in air pollution-health associations between 
the two cities. 

Despite this study's limitations, the potential to use spatial 
models of airborne pollutant concentrations to estimate exposure 
associated with asthma morbidity in Detroit and Windsor is 
encouraging. As new quantitative methods are developed to 
aggregate and incorporate health and environmental information 
collected at different geospatial scales, 57,58 they will have an 
increasingly important role in modeling, exposure estimation, and 
hypothesis testing in epidemiological studies. 59 Indeed, recent 
investigations highlight the importance of geographic scale in 
detecting racial disparities in health outcomes. 60,61 For most of the 
analytes included in this study, observed concentration gradients 
change across distances smaller than the dimensions of the zip 
code regions or postal sortation areas used to aggregate asthma 
events and pollution concentrations. We postulate therefore 
that analysis at a finer spatial resolution (i.e. census tract or 
block scale) can improve correlations between ambient air quality 
measurements and asthma health outcome variables, provided 
that a sufficient population resides within each geographic 
subdivision to maintain statistical power. 

In summary, concentration maps for the suite of measured 
contaminants demonstrate neighborhood scale spatial variability 
of air toxics distributions in Detroit and Windsor. Statistically 
significant positive associations between acute asthma events and 
N0 2 , VOCs, BTEX, and PM 10 concentrations were found at postal 
code scale resolution in one or both cities. Further study is needed 
to examine whether the air-health effects observed in 2008 
remain consistent in 2009 or whether the average value of the two 
sampling campaigns is a better proxy of exposure. Additional 
study is also required to investigate the degree to which the 
observed correlations were affected by the chosen level of spatial 
resolution, individual measurement points, synergistic effects of 
exposure to multiple air pollutants, community-level character- 
istics that can modify associations between air pollution and some 
health outcomes, or by other as yet unrecognized factors. 
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